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While the expansion rate ol a homogeneous isotropic Universe is simply proportional to the square- 
root of the energy density, the expansion rate of an inhomogeneous Universe also depends on the 
y—i , nature of the density inhomogeneities. In this paper we calculate to second order in perturbation 

variables the expansion rate of an inhomogeneous Universe and demonstrate corrections to the 
' evolution of the expansion rate. While we find that the mean correction is small, the variance of the 

correction on the scale of the Hubble radius is sensitive to the physical significance of the unknown 
■ spectrum of density perturbations beyond the Hubble radius. 
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t^J- , I. INTRODUCTION 

o ; 

There is no more fundamental physical quantity in cosmology than the expansion rate of the Universe. In recent 
ij years the present value of the expansion rate, Hubble's constant, has been measured with increasing accuracy 0]. 
JL i With the exploration of the Universe at rcdshifts of order unity, we now have information about the time evolution 
of the expansion rate Q- A most surprising result is that the time evolution of the expansion rate does not seem to 
be described by a matter-dominated Friedmann-Lemaitre-Robertson- Walker (FLRW) cosmological model. The usual 
explanation for the discrepancy is that there is a new component of the energy density of the Universe, known as 
dark energy, that determines the recent evolution of the expansion rate. Of course all indications for dark energy are 
indirect; they all involve some form of the time evolution of the expansion rate. 

Since the expansion rate of the Universe is of such fundamental importance, we must understand any possible 
effects that would result in an expansion rate different from the FLRW prediction. In this paper we study the change 
in the expansion rate due to perturbations of a homogeneous, isotropic, FLRW model. In particular, we perform a 
second-order calculation of the effect of inhomogeneities on the expansion rate. We only consider modifications to 
the expansion rate of a matter-dominated universe, although our results can be extended to a universe containing a 
mixture of matter and a cosmological constant. We find that the mean corrections are a few parts in 10 5 . 

The expansion rate of a perturbed FLRW cosmology has been discussed in many works (although we believe ours 
is the first complete second-order calculation). Hui and Seljak estimated the order of magnitude of the effect by 
considering a representative second-order term. They arrived at the correct mean order of magnitude of the result. 
Some recent works |5( suggested that small-scale contributions could give a large correction, producing an apparent 
accelerated expansion of the universe. In particular Rasanen p| suggested that due to ultraviolet sensitivity, one 
of the second-order terms could give large corrections depending on boundary conditions. As we will discuss in the 
conclusions, we find that if one employs the correct averaging procedure the result (at least for the terms we could 
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compute exactly using the second-order formalism) is well behaved in the ultraviolet and the term identified by 
Rasanen does not result in a large correction. 

In addition to calculating the mean value of the expansion rate, we calculate the variance about the mean value. 
In calculating the variance we uncover an interesting infrared effect. At second order the expansion rate has a term 
proportional to ipS/ 2 Lp, where y(x) is the peculiar gravitational potential, related to the density perturbations through 
the cosmological Poisson equation. If inflation is the origin of perturbations, then tp should be a Gaussian random 
variable with zero mean. However, if the perturbation spectrum on super-Hubble-radius scales is no bluer than a 
Harrison-Zel'dovich spectrum, then the variance of ip formally has an infrared singularity. 1 Now if the value of ip in our 
local Hubble volume is sufficiently large to modify the expansion rate, our perturbative expansion fails. Nevertheless, 
our results suggest that if the super-Hubble modes of ip have physical significance, then a nonperturbative extension 
of our calculation could yield a most important modification to the Friedmann equation. One might even speculate 
that a complete treatment of the effect could explain the observed time dependence of the expansion rate on its own 
and obviate the need for the dark-energy assumption. Since that conjecture is beyond the perturbative calculation of 
this paper, we will postpone discussion of this point to a subsequent communication Q. 

In the next section we discuss the general perturbative expansion. In Section IIIII we calculate the corrections to 
the expansion rate in terms of metric fluctuations, and then express the results in terms of the density perturbation 
spectrum. In Section llVI we present the numerical results. We then conclude, followed by two technical appendices. 

II. THE GENERAL FORMALISM 

In this section we describe how to treat the average properties of a perturbed Universe up to second order in the 
metric variables, including the effect of the inhomogeneous gravitational field on the homogeneous background field. 
We then discuss the proper definition of the average expansion rate and express the effect of inhomogeneities to second 
order as averages of density perturbations. 

By "average," we mean the average over a spatial hypersurface at a given time. Clearly an average defined in this 
way depends on the chosen coordinate system, i.e., the gauge. 

We will consider a Universe filled only by irrotational dust and choose the coordinates of an observer at rest with the 
dust (i.e., comoving coordinates), and with the same time coordinate for every point of the spatial hypersurface (i.e., 
synchronous coordinates) . This system of coordinates can be chosen if the Universe is filled by a single pressureless 
component. Since the pressure vanishes, the only force acting on the particles is gravity, and the comoving world lines 
coincide with geodesies. (See Ref. for a full discussion of this point.) 

We will call r the (conformal) time in this gauge, and x l the spatial coordinates, so that the metric has the form 

ds 2 = a 2 (r) [-dr 2 + 7y (r, x l )dx l dx j ] (1) 
We will perform the averages on constant-r spatial hypersurfaces. 

A. First order 

The goal of this paper is a second-order calculation, but before embarking, let us recall the familiar linearized 
first-order result. The energy-momentum tensor and metric are expanded to first-order as 

V = TW+rW (2a) 

= g$+g$, (2b) 

where the superscript (r) denotes the r-th order perturbation. By definition T$) and gj$ are homogeneous and 
isotropic. 

In the synchronous gauge, to first order the metric may be written in terms of a set of perturbation variables 
consisting of two scalars (ip^, X^)j a vector (x^); and a tensor (Xy*)- The vector xj 1 ^ is transverse (d t x\ 1 ^ = 0) 

and the tensor x\] is symmetric, transverse, and traceless (x\j — Xj\\ <9 l X« = 0, X i = 0)- in terms of these 
variables the metric is 

Hi = (l " 2^ (1) ) Sij + A,'X (1) + d iX f + djxP + xg\ (3) 



1 Presumably this infrared singularity is cut off because inflation did not last for an infinite period of time. 
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where Dy = didj — gV 2 5y. 

The metric perturbations tp^ and x may be expressed in terms of the peculiar gravitational potential <£>(x), 
which is related to the first-order density perturbation, by the Poisson equation, 

VV=yaV 0) <5 (1) . (4) 

Ignoring metric perturbations that decay in time, for a matter-dominated Universe jM 1 ) and x are given by (see 
Ref. 0) 2 

V^W) = ^(x) + ^W(x) (5a) 

X (1) (x,r) = -^rV(x). (5b) 

To first order we will neglect vector modes, as they do not arise in conventional perturbation-generation mechanisms 
such as inflation. We will also assume the tensor mode amplitude is small. 

In calculating spatial averages we will also require ^Fy, where is the determinant of the spatial metric. To first 
order 

V7 = 1 + ^7 (1) (x, r) = 1 - 3V (1) (x, r) = 1 — 5p(x) - yVV(x). (6) 

Whether calculating to first order or second order, when calculating the spatial average (• • • ) of a quantity 0(t, x l ) 
one must fix a system of coordinates in which to express 0(r, x 1 ), and then integrate with the proper integration 
measure over d 3 x at fixed r. For all orders we adopt the definition 3 

(0)(T)s /^V7^£fcg!) t (7) 

J d 3 x \J j(t, x l ) 

where the domain of integration is some large volume. If is a homogeneous quantity {i.e., it does not depend on 
x l , but only on r), then simply (0(°') = O^ \ In the first-order calculation, if is already a first order quantity 
we may set 7 = 1, and (O^) becomes 

To first order, the Einstein equations are 

G$+G#=^(t£>+T&>), (9) 
which yield upon averaging the 00 component 

G^=k 2 {T 00 }-(g^). (10) 

Here we see that kT 2 (Gqq \ may be interpreted as an extra component to the stress-energy tensor. 4 

For the unperturbed FLRW cosmology, the expansion rate a' /a 2 is found from the — component of the Einstein 
equations. For the first-order perturbed FLRW model 



2 



3(^1 =n 2 {p)-a- 2 (G^), (II) 



2 The perturbation variable <f> in Ref. 8| corresponds to our variable -0, and vice versa. 

3 Note that this definition is appropriate for scalar quantities. It will also be used to average Goo which is a scalar under the residual 
gauge freedom. 

4 This effect is sometimes incorrectly described as a "backreaction." Technically, it is not a backreaction; while the inhomogeneities 
modify the expansion rate of the Universe, they are not produced by the expansion of the universe. 
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where p = p( ' + pW = a 2 Tqq. This is the basic point. In a perturbed FLRW cosmology, a/a — a'/a 2 is not 
k 2 (p(°))/3. 5 However, one must be careful to find the true variables that describe the evolution of the averaged 

background. It is not clear that the quantity a'/a 2 describes the physical Hubble flow. The correct quantity to 

describe the Hubble flow may be found from the evolution of (p) . 

We know that for a homogeneous isotropic Universe there are two independent equations that govern the dynamics 

of the expansion. So we have to find the two independent equations that describe the evolution of averaged physical 

quantities such as like (p). For the unperturbed model we may augment the — equation with the continuity 

equation, D^T^ = 0. This gives for a perfect pressureless fluid (dust) 

-p'=p = -6p, (12) 
a 

where 9 = D^u^ and is the fluid four velocity. Eq. (|12|l is true to all orders. For an unperturbed FLRW Universe 
with scale factor a, this immediately gives p' / p = —3a' /a, which results in p oc a~ 3 . 

For the perturbed model, we want to expand Eq. I|12fl to first order and then average it. This will give us the 
effective scale factor for the dilution of the matter density. In addition to the expansion of p, we expand 9 as 

e = e^ + e {1 \ (13) 

and obtain 

~ + P V>') = ( (V°> + 0«) (>) + p«) ) . (14) 

Using the fact that (p^') = (p' 1 ')' we are left with 

^ = "<*>• (15) 
a (p) 

This tells us that in a first-order perturbed Universe the average matter density is diluted with expansion according 
to (p)~ 1 d(p) / dt = —(9), which in general is not equivalent to (p)~ 1 d(p) /dt = — 3a/a. 

An alternative way to express the result is to define a new scale factor, ay, such that (p) oc ay 3 . It is clear that av 
must be defined by 

? - >• (16) 

And so the true scale factor is ay and not a, at least for the dilution of matter (see, e.g., 0). 
The physical quantity of interest is (9) . It is given by 



(9) =0(°) + (0( 1 >) =3— + (9W). (17) 
a z 

We will define a 69, and express 9 as 



Two important notational points: We have defined H = yj k 2 (p)/3 (not as a'/a 2 ); and we have defined {89) as the 
difference between (9) and 3H (not as the difference between (9) and 9^°'). 
To calculate (69) /?>H, we can express Eq. (|llfl in the form 



(p) \^oo / \^00 / 

5 3^H ^ H 1_ 6^l' (19) 



5 Here the overdot stands for d/dt, where t is related to r by adr = dt. 
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where of course the second equality holds if the corrections are small. Combining Eqs. i|18|) and (|19|l . we find 

(59) (6>«) ( G oo 



(20) 



3H 3H 6a 2 H 2 ' 

It will turn out that in the first-order calculation gffi will appear only as a spatial gradient in the final expression 
for (59). This means that the physical result is insensitive to the choice of the normalization of (g t t))i since we could 
always add a constant to g$ to make (g$) anything we please. The cancellation of the non-gradient first-order 
terms in Eq. H20JI seems accidental, but in Appendix E] we derive an expression for (59) that at first order only 

contains spatial derivatives of g$ , which do not change if we shift gfy by a constant (or even a time-dependent 
renormalization) . 

B. Second order 

The spatial metric 7^ is expanded up to second order (neglecting first order vector and tensor perturbations) as 6 

m = (i - 2<A (1) - <A (2) ) 4, + Aj (x (1) + \x (2) ) + \ (9;xf + d jX ™ + x g)) . (21) 

The functions ip^ r \ X^ r \ X^[\ an d represent the rth-order perturbation of the metric. Vector and tensor modes 
have been included here at second order as they are dynamically generated by the non-linear evolution of purely scalar 
perturbations [8j. 

It will turn out that ^>' 2 )(x, r) is the only second-order term for which we will require the explicit form. For a 
matter-dominated Universe the metric perturbation (x, r) is obtained similarly to Ref. @ 

V/ 2 W) = -y V % ~ + ^ ((VV) 2 ~ • (22) 

The first term arises from a primordial epoch of inflation and can be computed as follows. There exists a second-order 
extension of the well-known gauge-invariant variable £, the curvature perturbation on uniform density hypersurfaces. 
In terms of H = a' /a, to first order it is given by C (1) = + HS^p/p'. To second order, Q = C (1) + (1/2)C (2) , 
where C,^ remains constant on superhorizon scales in the case in which only adiabatic perturbations are present. In 
standard single-field inflation, C,^ is generated during inflation and its value is given by ~ —2 (C^) [Tol Hl|. 
Since in the synchronous gauge and on superhorizon scales — — 5^/3 and C}- 2 ^ ~ tp^ 2 \ one readily concludes 
that the primordial contribution to the second-order metric perturbation ij)^ is given by the constant term — 50</? 2 /9 
and it also propagates to the coefficients of the term proportional to r 2 . 

In the second-order calculation care must be taken in defining the spatial average. Again, if O is a homogeneous 
quantity, then simply (O^) — O^K HO is already a second-order quantity, then we can take 7 = 1, both in the 
numerator and in the denominator. But in the second-order calculation of (O^) we must remember to include s Fy. 
Using Eq. ©, 

(1) /rf^y^^OWfoa:') = Jd 3 x (1-3^)QW(t,x*) 
{ ' ftPxy/^r, x 1 ) fd*x (1-3^(1)) 

where we have introduced the notation (• • • )i for a first-order term to denote the spatial average with the factor 
V7 = l, i.e., (G (1) )i(r) = fcPtcOWfax^/ft&B. 

So the correct second-order averaging procedure gives 

(OW) = O (0) (24a) 

(©(!)) = (0( 1 )) 1 -3(^ (1) (1) ) +3(?/> (1) )i(0 (1) )i (24b) 

m (24c) 



From now on spatial indices will be raised and lowered by the background metric Si 



6 



As an illustration of the subtleties in averaging, if at first order (gjiv}i = 0, where the first-order averaging prescription 
is defined in Eq. JSJ), at second order the averaging prescription for a first-order quantity is defined by Eq. Q24bfl . and 

(gfy) need not vanish. 

Now we perform a second-order expansion. We first expand the energy-momentum tensor and the metric as 

V = TW+rW + rcg (25a) 

9v> = 9$ +9$ +9$, (25b) 

where again the superscript (r) represents the rth-order perturbation. Again we will take TjxD and g ®) to be ho- 
mogeneous and isotropic, so T$) — (T^)) and g$ — (g$)- Note that this is not equivalent to the statement that 

T$ = (Tftu) and g$ = (g^ v ). 

Expanding the Einstein tensor G^„ to second order yields several types of terms. The first type are the zeroth-order 

terms, denoted Gj2, which only involve gffl . G^} will be homogeneous since gj$ is homogeneous by definition. Then 

there are first-order terms, G^J, which involve a single power of g$, possibly combined with g$- There are two types 

of second-order terms. The first type, denoted by G^\ involves squares of g$ ■ The second type of second-order 

terms, denoted G$, involves gffl- 

Now we consider the Einstein equations 

Gffl + + Gg? + G$ = - 2 (rjS? + T® + T$) . (26) 
We proceed by averaging Eq. (JSDJ. This yields 

G oo = ( T oo) ~ (poo + Goo^ + g oq ) • ( 27 ) 

From Eq. H27[) we see that now kT 2 (Gqq + Gqq 1 -* + Gg 2 ^) may be interpreted as an extra component to the stress- 
energy tensor. To second order, the — component of the perturbed FLRW model gives 



a 



2 



3(^1 = « 2 (p) - a- 2 GW + G&» + G$ ) . (28) 



Again we see that a/a = a' /a 2 is not n 2 p(°'/3. 

Now let us consider the second-order expression for the evolution of (p). We again use the continuity equation, Eq. 
(|T2*|l . and expand p and as 



p = p^+p^+pW 



e = 0(°) + o(i) + o(n) +e (2)_ (29) 
The result is a bit more complicated than the first-order result: 

I ( p (o)/ + p (u> + p (2)/^ = _ ( ^(o) + 6 (d + 6 (n) + 9 wj ( p «>) + p( i) + p wj ) m 

Using the fact that (pW') = (pW)' - 3(/jW)i(^W')i + p {1) ) , the left-hand side of Eq. © becomes 

- + P W + P M>) = -(P)' - -(^ (1) ')i(P (1) >i + -(^ (1) 'P (1) >. (31) 

a \ I a a a 

We may also express the right-hand side of Eq. (|30|) as 

_ ((0(0) +0(1) +0(H) +0(2)) ( p (0) + p (l) + p (2)^ = + (0(1)),^!))! _ (p«0«). (32) 

Equating Eqs. (|31(l and Ij32(l we obtain 

i(p>' - -{^'hip^h + = + (tfW)!^ 1 ))! - (P (1) ^ (1) >. (33) 
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Since to first order a simple calculation yields = — 3-0' 1 ) /a, we are left with 



a (p) 



7^ = -<*>• 



(34) 



This tells us that as in a first-order perturbed universe, in a second-order perturbed Universe the average matter density 
is diluted with expansion according to (p)~ 1 d(p) jdt = —(0), which in general is not equivalent to (p)~ 1 d(p)/dt = 
—3d/ a. Note that it was crucial to define the averages with the factor; otherwise we would not have discovered 
the right quantity to describe the scaling of (p) . 

Again, the physical quantity of interest is (9), given by 



(9) = 0<°> + (0* 1 )) + (9^) + (9^) = 3^ + + (0( n >) + (9^). 



(35) 



We emphasize again that H = \J n 2 (p) /3 (and not a' /a 2 ). 

The goal of this paper is to calculate (89) /3H. To do so, we can express Eq. (12811 in the form 



n 2 (p) 




G, 



(i) 



G, 



(ii) 



G 



(2) 



G, 



(i) 



6a 2 H 2 



72a 4 H 4 



where of course the second equality holds if the corrections are small. Using Eq. H3(jfl . we find 



3H 



(ew + 9^) + e<») (4o + 4o 1} + 



G 



3H 



6a 2 H 2 



72a 4 ff 4 



(36) 



(37) 



(2) 

It will turn out that g^J will appear only as a spatial gradient in the final expression for (89). This means that 

(2) 

now the physical result is insensitive to the choice of the normalization of (g^J), since we could always add a term 
that is spatially constant. The cancellation of the non-gradient second-order terms in Eq. I|37|l seems accidental, but 
again, in Appendix [S] we derive an expression for (89) that explicitly only contains spatial derivatives of g$ , which 

(2) 

do not change if we shift g^J by a term that is spatially constant (or even a time-dependent normalization) . 



III. COMPUTATION OF (S0)/3H IN THE SYNCHRONOUS GAUGE 

The relevant quantity which we want to calculate is (89) /3H. The simplest way to proceed is to perform the 
computation directly in the synchronous gauge. The synchronous coordinates are more physical for our purposes. 
If the calculation is performed in the Poisson gauge it is necessary to perform a complex calculation (as we do in 
Appendix [B| to express the result in synchronous coordinates (see also Refs. p| and 0])- 

Moreover, in general for second-order calculations, the synchronous gauge turns out to be very convenient, since 
here the scalar perturbations do not have nonlocal terms that appear in other gauges such as the Poisson gauge. Also, 
for matter it M in this gauge takes the trivial form u M = a _1 (T)(l, 0), so the only terms that have to be computed are 
the Christoffcl symbols in the covariant derivative 

9 = = d»u» + r>* = 3^ + 5T ° oa + 6r °\ (38) 

p a 1 a 

Moreover, in the synchronous gauge, to second order STq° = 11]. Finally, the computation of the perturbation of 
9 with respect to its background value consists of finding just the trace of rj, ■. From [D| , 

3H ~ 3H {d 0i) ~ 3Ha V M + 6 °7 

= ^ - - 2^ (1 ¥ 1} ' + ^VV 1} V 2 x (1) ' - \x {X)M X [1) ', k ) 

= -rU /--VV - — ^VV + — <P " — (VV) 2 - — V'^V a) • (39) 
a 2 H 2 \ 9 v 2T * 2T^ 63 v ^' 189^ / v ; 
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In the last equation we expressed the result in terms of the initial perturbation tp, using the time evolution of the 
relevant perturbation variables in the synchronous gauge given in Eqs. JSJ and 1221) . 



The simplest way to find Gqo is to realize that a 



^oo 



(ii) 
oo 



G 



(p 



(i) 



P 



(2) 



and use the explicit 



solutions for and from Ref. but using for ip^ the expression augmented with the tp 2 as in Eq. (|22() . In 
the matter-dominated universe, 



G, 



(i) 



G, 



(ii) 



G, 



(2) 



6a 2 H 2 



a 2 H 2 



-vV 



10 ~ 55 



5rf 
126 



G, 



1 



72a A H 4 



a 2 H 2 72 



(vV)(vV>- 



(40) 



The value of 59 /3H is found by summing Eqs. (EH and (|4*0"|) . 

We must now perform the appropriate spatial average. In 86 there is only one first-order term; the rest of the terms 
are second order. The first-order term must be averaged using the procedure of Eq. I|23ll . while the other terms are 
averaged using y^y = 1. The result is 



3H 



a 2 H 2 



3 

23t 2 
216 



(41) 



A. Evaluation of (89) in terms of the matter power spectrum 



We now proceed to express the averages in terms of the matter power spectrum. The procedure is to fix a spherical 
domain of radius R with volume V(R). From Eq. (0) and the definition of (• • all of the averages in Eq. i|41|) 
involve integrals of the form 



1 



V(R) 



(• • • ) d 3 x. (42) 

V(R) 



For calculational convenience we will employ a Gaussian window function and assume V(R) is a spherically symmetric 
volume with volume clement dV = 4nr 2 exp(— r 2 /2R 2 )dr and volume V(R) = (2tt) 3 ' 2 R 3 . The Fourier transform of 
the window functions is 



W{kR) = V-\R) J d 6 x e' r exp(zfc • x) = e'" " I*. (43) 

Of course as kR -> 0, W{kR) -> 1. 

We wish to evaluate the typical expected value of 9 averaged over this sphere. By "typical expected value" we 
mean the ensemble average. The metric fluctuation p is treated as a Gaussian variable with zero mean (of which we 
know the iV-point correlation functions) that takes random values over different "realizations" of volumes V(R). In 
other words, we calculate the typical value of a quantity for a region of radius R as the statistical mean over many 
different similar regions. We will indicate this statistical average with a bar: (•••). 

We will express ip and its derivatives in terms of a Fourier integral, so 

d k rit.^ f d 3 k . ,9 f d 3 k 



v = J Wf n elk ' S ' = J Wf lkl n vV = " J Wf fc2 Vi etc - (44) 

The Fourier components ip? satisfy 

n = (45a) 

p^Tp~l = (2n) 3 S^(k 1 + k 2 ) P v (h) (45b) 
Wp?; = (27rf{8^\k 1 +k 2 )8^(k 3 + k A )P v (k 1 )P v (k 3 ) 

8^{h + k 3 )S (3 \k 2 + k 4 ) + S^ih + k 4 )S^{k 2 + fc 3 )] P v (fci)P v (fc2)} , (45c) 
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where P v (k) = \ft\ ■ From Eq. Q we can express P v (k) in terms of the matter power spectrum as 

9tt 2 



(46) 



where A 2 (/c,a) is the (dimensionless) power spectrum of the matter density fluctuations. 

Let us first consider (<p) 1 and (V 2 </j) r Clearly from Eq. (|45a|l . (ip) 1 = and (V 2 ^^ = 0. However this does not 
imply that ((p) 1 = or (\/ 2 ip) 1 — over any individual volume of radius R. The question is the magnitude of typical 
departures from the mean values, which corresponds to the statistical variance of our quantities. As we will show, it 
is intuitively clear that if the radius R is big enough, this variance will go to zero. So the effect of variance could be 
important. We will return to the calculation of the variances of the different terms after completing the calculation 
of the mean values. 



Next, consider ((p\7 2 ip). Passing to Fourier space, we have 

d 3 x d 3 k x d 3 k 2 f kf + k. 



V{R) V(R) (2tt)3 (2tt)3 



i(ki + k 2 ) • x 



d 3 ki d 3 k 2 fkf + k 



(2tt) 3 (2tt) 2 



mn 



TW 



( fci + k 2 R) 



Making use of Eq. 145b|) , we find 



d 3 k 



k 2 P v {k). 



(47) 



(48) 



Using Eq. (|46(l . we can express (ip'V 2 ip} in its final form. In the same manner we find the means for the other terms. 
The result is 



(<P)i = 



(V 2 ^>! 



\a 4 H 4 [ ^ 2 { k,a 



-a 4 H 4 



dk 
k 3 



A 2 (fc,a 



4a 4 £T 4 f°^A 2 (fc,a 
4 Jo k y 



(¥>}i<V 2 ^ 



-a 4 H 4 



a 4 H 4 



dk 



dk 

F 



A 2 (fc,a 
A 2 (fc,a 



l a i H i r^A 2 (k,a 
4 Jo k 



W 2 (kR) 
W 2 {kR). 



(49a) 
(49b) 

(49c) 
(49d) 
(49e) 
(49f) 
(49g) 
(49h) 



Of course in the limit R — + oo, Eqs. | |49g| | and lf49h|) vanish. Therefore, the entire second-order calculation gives 

where we have used r = 2/aH appropriate for a matter-dominated Universe. We will give numerical results in the 
next section. 

Now for the variances of selected terms. The variance is defined as 



Var 



[<•••>]=(<•••>-<•••>)' 



For instance, 



Var^J = \a 4 H 4 J^ ^A 2 (k,a) W 2 (kR) 
Var[<VV)!] = \a 4 H 4 J°° * A 2 (fc, a) tU 2 (fci?). 



(51) 

(52a) 
(52b) 
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The variances of other terms are more complicated, but straightforward to derive. For instance, 

-l 2 

W 2 (\h + k 2 \R) P v (Ai)P»(fe). (53) 



Var 



13 



i 



(27r) e 



d%d 3 k 2 - t 



1 



13 - 

k 2 + k 2 + — kt ■ k 2 



The angular integrals can be expressed in terms of a filter function, defined in general as 



J {l) {k ll k 2l R) = J <W W 2 (^Jk 2 + k 2 + 2k 1 k 2 fiR 



(54) 



For a Gaussian window function the filter expression can be expressed in terms of incomplete T functions as 
J^ik^k^R) = J 1 d/i (j, 1 e -(fef+ fe ?+ 2 P fe i fc 2)-R 2 

exp [- (fef + fcf) i? 2 ] [r(l + 1, - 2R 2 k 1 k 2 ) - r(I + 1, 2i? 2 fc 1 A; 2 )] 



2'+ 1 i? 2 ('+ 1 )fci fc 2 

We will make use of the fact that for k\R — > and /e 2 i? ~> 0, J^°'(ki, k 2 , R) — ► 2. 
Using the Gaussian filter, Eq. (|53|) becomes 



(55) 



Var 



13 



j^ 4 ) 2 f If ^ ' a) I " |f ^ (fc2 ' a) [ (| + | + 2 ) J(0) M > R) 



26 /jfei fe 2 



9 \k 2 ki 



169 



+ j~ J (1) (ki,k 2 ,R) + -^J {2 \k l7 k 2 ,R) 



81 



(56) 



It is interesting that the integral is not well behaved in the infrared. As discussed in the next section, for a Harrison- 
Zel'dovich spectrum A 2 (k) ex fc 4 , so the term proportional to has an infrared divergence: 



Var [(<pV 2 <p)] 



J JO K l JO K 2 \ K 2 K l/ 



(ki,k 2 ,R) 



(57) 



We will discuss the importance of this term below. 

Now we turn to the variances of (V 2 <p) 1 and (V 2 y) (V 2 ip). They are given by 



Var [fa), (VV)J 



9 a 4 H 4 



| ^A 2 ^,^ 2 ^) | ^A 2 (fc 2 ,a)tU 2 (fc 2 i?) 



A 2 (fc,a)PF 2 (fci?) 



Var [(VV> (VV>] = 2Qa 4 i/ 4 ) 2 (|^A 2 (fc, a )^ 



(JsR) 



Note that Eq. 



also has an infrared divergence 

9 , 



Var[(^) 1 (V 2 ^) 1 ] 



(58a) 
(58b) 

(59) 



Of course we will be interested in the variance of the total expression Eq. (|41|) . not the individual terms. Of 
particular interest is the cross term of the infrared-singular parts. This will be the only one for which we will include 
the cross terms. The infrared-singular pieces appear in Eq. 141|) proportional to (ip'V 2 (p} — 3(<^>)i(V 2 <p)i. The infrared 
part of the variance of this term is 



Var[(^VV)-3(v>)i(VV)i]c 



dk 



2 A2, 



dk\ 



A 2 (fc 2 ,a) / ^A 2 (fc l7 a) 



(60) 
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It is straightforward to obtain the variance for the second-order, four-derivative terms. (The sum of these terms 
has zero mean.) It is given by 



Var 



(vV) 2 ) - <w y ) 



o fc i 
+J^(ki,k 2 ,R) 



^A 2 (fc 1; a) r ^A 2 (fc 2 , a) \j^(h,k 2 , R) - 2J^(k u k 2 ,R) 



(61) 



There is another potential contribution that will result in contributions to the variance similar to the terms we have 
found. Suppose we expand 59 to third order in perturbation theory. We can express 59 in the general form 



(59) = (Atp + Btp 2 +C<p 3 ) , 



(62) 



where A, B, and C are operators which contain derivatives. Then the variance of 59 will contains terms like (A(p) 2 , 
(Bip 2 ) 2 , and (Atp)(C(p 3 ). The first term is the usual cosmic variance term; it is well behaved in the infrared. The 
second term is singular in the infrared; it is given in Eq. (|60ll . The third term will be present, in principle it also 
will have an infrared singular part, and its value requires a relativistic third-order perturbation calculation. However 
there is no reason for the infrared singular part of the AC term to cancel exactly the infrared singular part of the B 2 
term. 

Finally we remark on the ultraviolet behavior of the corrections to the expansion rate. The second-order result 
for the mean, Eq. I|50|l . should be well behaved in the ultraviolet. As discussed in the next section, in the linear 
regime A 2 (fc,a) increases logarithmically with k for a Harrison-Zel'dovich spectrum, so the first term should be well 
behaved in the ultraviolet. The ultraviolet behavior of the last two terms are regulated by the filter function W 2 (kR). 
The contributions to the variance of the terms we have calculated, Eqs. (j52b(l . I|56|) . (|58|l . and (|61|l . all involve filter 
functions that regulate the ultraviolet behavior. However, we expect there to be terms that do not involve filter 
functions. For instance, if one performs the relativistic third-order calculation, one expects to find contributions to 
(69) from terms like ((V 2 </?) 3 ). The variance would then include evaluation of terms like (W 2 ip)((W 2 ip) 3 ) (an example 
of the aforementioned AC terms). These terms would include parts with a momentum integration unregulated by 
a filter function. If there are third-order terms with large numbers of derivatives bringing down large powers of 
momentum, then the variance might be sensitive to the ultraviolet behavior. 

In the next section we show the numerical results for the mean and variance of the corrections to the expansion 
rate. 



IV. NUMERICAL RESULTS 



In this section we present the numerical results for (69) /3H for a matter-dominated Universe. We will give the 
mean values, as well as the variances. 

For both the mean and the variances we express the power spectrum A 2 (fc,a) in terms of the transfer function 
T 2 (k). For a Harrison-Zel'dovich spectrum, the power spectrum is 



k x4 



A 2 (fc,a) = ^— j T 2 (k), (63) 

where A is the dimensionless amplitude, A = 1.9 x 10~ 5 . We will discuss the implications of other spectra. For 
our purposes the Bardeen, Bond, Kaiser, Szalay (BBKS) transfer function ^| will be adequate. The BBKS transfer 
function may be expressed in terms of a dimensionless parameter 

(64) 



r/iMpc" 1 ' 

where T is the shape factor, defined for a flat universe in terms of the baryon fraction 17 # and the total value of 
as r = ttoh exp(— — ^/2MIb/^o)- In terms of q, 

T(q) = hl(1 2 + 34 2 g 34g) [1 + 3.89g + (W.lq) 2 + (5.84g) 3 + (6.71,) 4 ] ~ 1/4 . (65) 

Of course at small q, T 2 (q) —* 1, while at large q, T 2 (q) — » q~ 4 ln 2 g. 

Also, in all expressions we make use of the fact that in a matter-dominated Universe H 2 (a) = H^a^/a? , where do 
is the present value of the scale factor. 
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k MAX (h Mpc" 1 ) 



FIG. 1: The present mean value, (59) /3H, in the matter-dominated Universe as a function of the ultraviolet cutoff &;max for 
various values of the shape factor V. 



Consider the mean, given by Eq. |(5 



m 

3H 



25 a , 

A 

24 a 



h Mpc" 1 \ p2 ' - 



d qq T\ q ) + 2 ^A 2 ( hMpc ~ X ) r 2 

4 a 



A" 

96 \a 



H 



Ha 



dq q T 2 (q) W 2 (rTq) 



dq q 3 T 2 (q) W 2 {rTq), 



(66) 



where r is the dimensionless size of the region, r = R/h~ 1 Mpc. We will present results for R = Hq 1 (so r — 3000) 
and the last two terms are negligible. While we have indicated the range of integration from q = to q — oo, in 
reality to employ Eq. H63(l there is a cutoff on the maximum value of the integral. The ultraviolet cutoff arises because 
density perturbations become nonlinear. The mean value is very insensitive to the ultraviolet cutoff. The integral for 
the mean value receives most of the contrib ution in the decade between 10" 1 < q < 1. 

The result for the present mean value of (86) /ZH is shown in Fig. ^for various choices of V. It scales as a/ao = 

Now consider the variance about the mean value. First consider the variance of the linear term, Eq. I|52b|) . This 
term contributes to the variance in (89) /3H an amount 



VVar [(86)} 



3H 



5 

9a 2 H 



WvaxKW)!] 



linear term 

Defining a dimensionless wavenumber x = qTr, the term becomes 

2 



dk 



A 2 (k) W 2 (kR) 



1/2 



(67) 



VVar [(86)} 



3H 



linear term 



5 (hMpc 

6 l H 



-l 



a 



(IQ 



dx x 3 T 2 (x/Tr) e~ x 



1/2 



a 100 

a r 2 



(rr>l), (68) 



where the last expression holds for xmax = ^max-R ^ 1- The result is shown in Fig. [21 Because of the window function, 
the results do not depend on /cmax, the ultraviolet cutoff (so long as it is greater than about k = O.lh Mpc -1 ). The 
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R (h-« Mpc) 



FIG. 2: The present contribution to y'Var [{S6)]/3H from the indicated terms in a matter-dominated Universe as a function 
of R, the volume sampled. The results are for &max = 0.1 and T = 0.2. On the logarithmic scale, the dependence on V and 
fcmax is not noticeable. We emphasize the point made in the text that the calculation for the the second-order, two-derivative 
term has an infrared cutoff of the present Hubble radius. 



result is also well behaved in the infrared. Of course as R — > oo, the variance of the linear term approaches the mean 
of the linear term, which is zero. 

We can see comparing Fig. □ and Fig. that the variance in the linear term dominates the mean value out to 
distances of about 3 x lO 3 ^ 1 Mpc. 

Next, consider the contribution to the variance of a typical second-order, four-derivative term, Eq. This term 

contributes to the variance in (66) /3H an amount 



VVar [(89)] 



1 



3H 

second-order, fc 4 

Again, with the same dimensionless variables 

1 



27 a A H 4 



Var 



(69) 



y/Vax[(80)] 



3H 



second-order, k A 



= -A 2 
3 



h Mpc 



dxi / d,X2 x 3 , 
it Jo 



x T 2 {xi/Tr) T 2 {x 2 /Tr) e'^e^ / d^i e -**x*2i* (l - M 2 ) 



1/2 



(70) 



The results are also shown in Fig. This term scales as r~ 4 . The result is well behaved in the infrared and an 
infrared cutoff does not need to be introduced. 

Finally, consider the variance of a typical second-order two-derivative term, Eq. 1)56(1. As mentioned in the previous 
section this term is not well behaved in the infrared. We first calculate the variance by introducing an infrared cutoff 
which we will take to be the wavenumber of the present Hubble radius. The second-order k 2 term contributes to the 
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variance in (59) /3H an amount 



VVar [(69)} 



3H 



25 1 
2Ja 2 H 2 



<Var 



13 

((p\7 2 ip) + — ((fjip'i) 



(71) 



second-order, k 2 

In terms of the dimensionless variables x and r, this term contributes to the variance an amount 

2 



VVar[(50)] 



3H 



second-order, fc 2 



24^ 



ft Mpc" 



IMAS 



IMAX 



dx 2 T 2 (xi/I>) T 2 {x 2 /Tr) e^e^ 



2G 



2G 
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dfi e' 2xiX2t " M- + — + 2xia; 2 + — + — x\\x + —xiX 2 fJ,' 
Vaf2 9 9 81 



1/2 



(72) 



As a result of the infrared behavior, for this term we must introduce another parameter, xmin = ^min-R, as an infrared 
cutoff. Its value is shown in Fig. [21 as a function of R. As expected, it has a r~ 2 behavior. 

Now we turn to the infrared-singular part of the second-order two-derivative term. The contribution is 



VVar [(59)] 



3H 



second-order, k 2 



25 1 
27 



(73) 



We will evaluate the integral by defining some X* = such that xmin <C x„ < 1. We will be interested in large 
values of r, so if x <C 1, then ir/IY <C 1, and T 2 (x/Tr) can be set to unity. Using Eq. I|6U|) . this term may be written 

as 



3i? 



second-order, fc 2 



25 ,2 

12 



ft Mpc' 



«o 



iz^ In ■ 



-,1/2 



'MIN 



(74) 



If we take fc* = fc# where kn = H = 3000 x h Mpc 1 is the Hubble wavenumber, then evaluating the term for r 
corresponding to the Hubble radius we find 



VVar [(59)} 



3H 



second-order, k 2 



25 j9 a 

— A A 

24 a 



In- 



>-MIN 



1/2 



25 . a 



3d ao 



(75) 



We started with a perturbative expansion, and for that expansion to be valid requires ip < 1, so our perturbative 
analysis will break down for \/Var[(<5#)]/3if ~ 10~ 5 . 

Of course the total variance of (59) /3H includes cross terms from the various contributions. However we have seen 
that the linear term will dominate. 



V. CONCLUSIONS 



Our results are illustrated in Figs and The mean corrections to the expansion rate are a few parts in 10 5 . If 
the correct prescription for dealing with the infrared singular nature of Var[y>] is to employ a cutoff of the order of the 
Hubble radius, then the variance about the mean is small, and dominated by the first-order term. If, however, the 
super-Hubble modes are physical, then the variance is dominated by the infrared singular part of the second-order 
corrections. Restricting ourselves to Var[y>] < 1 where the perturbative calculation is valid, then the variance will be 
small, of order 10~ 5 . 

In Ref. |5j , Rasanen considered what amounts to the variance associated with the second-order term proportional to 
W 2 ipW 2 (p and suggested that due to unknown boundary conditions it may lead to a large contribution to the expansion 
rate. We claim that when properly averaged, this term is subdominant. 

In closing, let us speculate that corrections to the expansion rate will include a term like ipS/ 2 ip even when ip 3> 1. 
Returning to the infrared-singular part of the second-order variance, Eq. (|75|) . 



= 3.76x l - 10 -ln-^-. (76) 
an Kmtn 



VVar[W] 
3H 
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For the variance to be of order unity, the perturbation spectrum would have to extend to a factor of exp(6 x 10 18 ) 
(10 18 - 8 e-folds!) times the present Hubble radius. 

Rather than a Harrison-Zel'dovich spectrum, if we assume a slightly red spectrum so that A 2 (fc) cx k 4 ~ e with 
< e -C l, 7 then the logarithmic term in Eq. (|76[) is replaced by e _1 (fc^r/A: m i n ) e . Now this will give unit variance if 
In fcff/fcMiN = (43.3 + lne)/e. For instance, if n — 0.94 on super-Hubble-radius scales, then a variance of order unity 
is obtained if the perturbation spectrum extends 676 e-folds beyond the Hubble radius. Since the present Hubble 
radius corresponds to a scale that crossed the Hubble radius about 60 e-folds before the end of inflation, if inflation 
lasted more than 736 e-folds with a super-Hubble-radius spectral index of n — 0.94, then the effect of super-Hubble- 
radius perturbations would have a large effect on the expansion rate of our Hubble volume. Speculation about the 
significance of this result will follow in a separate paper [6j . 
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APPENDIX A: ON THE PRESCRIPTION FOR T$ AND g$ 

(2) 

The aim of this appendix is to demonstrate that (9) does not depend on the value of (g^J) since second-order 
irms will only enter as sj 
constant-r hypersurfaces: 



terms will only enter as spatial derivatives. For the demonstration, we first introduce i?*-, the extrinsic curvature of 



if 



1 



1 2 

In either the synchronous or comoving gauge one can see that (see Ref. @) 

0(D + 0(ll) + 0(2) = ^ = t 

a a 

One can now write the G° term in a convenient form (see Eq. (4.3) of 



G, 



+ C$> + = \ {lU + 3 K + , 



where 3 1Z is the intrinsic curvature of the three space with metric 7^ . At second order 3 1Z is 



j ,lk 



+7 ,i (7 3 ,k ~ 7jfc 
Using Eq. l|A"2)l . Eq. also can be written as 

1 



r-(2) _ ■ 



4a' (d<n + + + 3 + 3 k^ + 3 k^ + (^(D) 2 - 

Recall now that in any of the definitions adopted, Eq. i|37[) is always valid. So independent of (gffi), 



(9) = 3H - 



(1) _L/o(ll) _l G (2) 



Goo + ^00 
2a 2 H 



(10 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



In the the usual parlance, this corresponds to a scalar spectral index of 1 — n = e. 
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Using Eq. (|A5(1 . we obtain 

(0) =3 H- ^ [( 3 ^ (1) + 3 ^ (11) ) + ( 3 ^ (2) ) + ((^ (1) ) 2 - • ( A7 ) 

We note now that the only piece that depends on gffi is the one contained in 3 TZ^ 2 K And, as we can see from its 
explicit expression 

3^) = l( 7 (^ fe _ 7 (^) j (A8) 
(2) 

this is a spatial gradient, so independent on (g}iJ). From this formula it is also clear that 6 has no term with 
zero powers of k, since the spatial curvature 3 1Z always involves spatial derivatives (and the same is true for the 

- terms). 

If one computes (9) using Eq. (|AT|) . one obtains the same result as Eq. I|41|l . thus our results are independent of 
(gffi)- By the same argument, in a first-order calculation we could add a constant to gftj. 

APPENDIX B: 9 IN THE POISSON GAUGE 

In this appendix we repeat the calculation in the Poisson gauge in order to check our results and compare with other 
results in the literature for 8 in the that gauge. 8 The Poisson gauge is a second-order generalization of the longitudinal 
(or Newtonian) gauge. Since this is a very involved way to do the computation, and since it is only intended to be 
a check, we do it only in a partial way by just computing the corrections + O^- 11 ^ + 6^ [see Eqs. (|18|) and (|39l) ]. 
Furthermore we compute exactly the terms with two spatial derivatives, while as for the four derivatives terms, we 
only check that they have zero statistical mean (that is, we disregard total spatial gradients). 

The metric is given at second order in this gauge by 9 



ds 2 = a 2 P ( V ) {- (l + 24> (1) + 4> (2) ) dr, 2 + Jpdy'dq + [(l - 2^> - ^) + \ 



dy l dyi . (Bl) 



(2) 

2 Ay 

(2) 

Here xlj is a pure symmetric, transverse, traceless tensor degree of freedom. We will denote the derivative with 
respect to conformal time rj as d v . We will also introduce here a cosmic time z, linked to r\ through the relation 
dz = a dr/, and we will denote the derivative with respect to z with a d z symbol. The spatial derivatives with 
respect to x l will be denoted as before, with a simple d % or ,l symbol. The derivatives with respect to y % instead will 
be explicitly written as d/dy l or with the subscript d\ y y Recall also, as well known, that the i — i component of 

Einstein's equations (in the absence of anisotropic stress) imposes (f>^ = ip^- 1 ' , and moreover, the evolution equations 
in matter domination give 

4,^(z,y i ) = ^ 1 \z,y i ) = i P (y i ). (B2) 
The four velocity of the fluid here has a more involved form 



it 



- h i 1 -* + ^ 2 + <4 9(y)i ^ _ ^) ' (B3a) 

u% = ^ P {-m- p d ^ + —)> (B3b) 

where Hp = Op^d^ap. Here v^ 1 is the second-order contribution to the three velocity. We will need only its 
divergence (which is obtained taking the divergence of the — i Einstein equation): 



2 / , , v 2 ,^ 2 )'\ 

' V? y) *< 2 > + -^g— . (B4) 



3H P \ {y)V H 



Actually, the other results we are aware of in the literature use only the first-order metric. 

(2) (2) 

Also here, for our purposes the vector and the tensor x\j wm never enter in the computations. 
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Now, taking these expressions and taking Christoffel symbols from Ref. > we must compute 

9 = 9X + T« a u» = d ri u° + d iy) iu* + T^u". (B5) 

This gives us 

Here Hp is defined as Hp = dp d z ap — dp d^ap = dp Tip. The piece coming from c*( y )jW l in Eq. (|B5(1 contains first 
order (<9( y ) ^u^ 1 ' J ) and second order {d( y )iU^ 1 ) terms. The first order term [see Eq. l|B3b|) ] is equal to 

^ (1)l = -^k v ^ (B7) 

and it will be important in the next steps of the computation, since it will produce second-order quantities. The term 
d{y)iU^ 1 [Eq. (|B4|) ] is a spatial gradient of second-order quantities. For the purposes of this appendix we have to 
keep the first two pieces in Eq. (|B4|I . and moreover we have to consider also the two intrinsically second-order terms 
<fP^ and ip( 2 \ Note also that, in their time evolution they contain non-local quantities (&o and \&o of Eqs. (6.8) of 
Ref. 8j), i.e., they are defined through 

v (y) e o = ^o-\d (y)lV d\ y) ^ 



(y)*° = ~\ 



( V (y)^) 2 - %b%)*Ky) a (y)H • (B8) 

So the spatial derivatives of these quantities will produce not only spatial gradients, but also terms who do not have 
zero statistical mean. In the end, we keep them in the calculation as 

where "• • •" indicates terms with four spatial gradients (that we will omit). 

So, the relevant expression for 9(z,y t ) is obtained by Eqs. I|B6|) . I|B7|) . and (|B9|I 

OM) = 3H P - 3H P , + "Jf + 2 °^g V) " - ^ + ^ 



3H<P^ 2V2 y)V , V 2 y) ^ 



3a 2 p H P 3a 2 p H P 3a 2 Hj 



(BIO) 



This is the result in the Poisson gauge in coordinates (z,y l ). Note that here there are terms which do not vanish 
in the infrared limit, proportional to ip and in ip 2 (they will disappear going back to the synchronous and comoving 
coordinates, as noted in Ref. |12|'). 

In order to compare with the results already obtained in the synchronous gauge, we have to express our quantity 
9{z,y l ) as 9(t,x l ) and average it over a volume in coordinates x % at constant r as in Section ITTll First, we have to 
change time from z to t in all the quantities. This is relevant only for the zeroth-order and first-order terms, and not 
for the ones which are already second order. So the quantities in which the time has to be changed are Hp, ap{z) 

and <p(z, y 1 ). The first two, cip(z) = (z/zq) 2 ^ 3 and Hp = 2/3z, have to be expanded up to second order in St, which is 
defined by t = z + St and expressed in terms of the scale factor and the Hubble rate defined by a comoving observer: 
a = (t/z ) 2 ^ 3 , H = 2/3i. Instead, for the variable (p{rj,y l ), ip(-q,y' 1 ) — (p{r,x l ) holds since d n (p(r],y l ) = 0. In order 
to find the change of the time coordinate from z to t one knows that u^d^ — dt, and that dtt = 1. In this way one 
obtains an iterative equation 

(l + + °) ± (z + St W + ftCO) + *d (y) i (z + St&) = 1, (Bll) 

where w M = ap 1 (<5o + v + v ^)- For zeroth, first, and second order, the equation becomes 

^ dz 

dz 

= v^ + ^St^ 
dz 

o = w (2)o + w (i)o * 5t m + * st m + v p)i d , Y)i 6tw, (B12) 

dz dz ^ ' 
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which lead to 

&t<V = (ft 



ftP) = -IvH+^gL + l fdt'^. (B13) 



2^ A5a 2 H 2 . 
Now using these equations, one can change the time in Eq. (|B10|) . obtaining 



11 8yMy (y )|i 3ii^ 2 > 3^ 2 > ; %H 2 f 4> {2) 
45a 2 H 



9(t, y*) = 3H+ \^Z {y) " - —7T V + — + d (y) y- ( B14 ) 



As said before, this transformation (i.e., passing to synchronous time t) has eliminated the two terms in <p and (p 2 
that do not vanish in the infrared. We wrote it in this form in order to compare with a previous result in literature 
(Eq. (10) of Ref. 0). As we are going to show, we agree with that result, but we added here the contribution of 
intrinsic 10 second-order terms (j)^ and ip^ 2 \ This is crucial if we want to recover the result obtained in Sec. IIIII in 
the other gauge. Here <9( y ) iU l contains the terms 

2V 2 sw 2wV? sw o / \ V 2 ,6^ V? ,ii (2) 

a ->-' —itfr-yf-ds K>*w + - w- - -£bt + ■ ■ ■ ■ (B15) 

Note that after changing time, a new term (which is not a spatial gradient) appears, namely —2ip'S/ 2 y ^ip/9a 2 H. 

Finally, one has to change also spatial variables from y % to x % . The change of coordinates is found by imposing that 
the metric g^ v in the coordinates (t,x l ) has goi = 0. Using the transformation given by Eq. 1|B13(1 . one finds that to 
first order the change in spatial coordinates is given by 

x*=y* + l-%£. (B16) 
3 a z H 

Since it is already first order, this change of coordinates is irrelevant for quantities that are already second-order. 
Moreover, the only first order term of Eq. (|B14|) is di y \ --i^ 1 - 1 ' = —2\7 2 y - ) (p/3a 2 H. Here, the change of spatial coordinates 

in the argument of <p(r, y l ) is not relevant for our purposes, since it produces only an extra term which is a four 
derivative spatial gradient of something. But passing from V 2 y ^ to V 2 produces a (relevant) extra term with four 
spatial derivatives, 

V J y) 9 = ~^VV ~ ! 4 Zi k - (B17) 



In all the other terms we can safely put d( y ) = d. After doing this, the result is 

2\7 2 ip 8tpX7 2 (p ^<fi' jk( P,jk 



0(t,x l ) = 3H 

3a 2 H 9a 2 H A5a 2 H 9 a 4 H 3 

3^ 2 > , 9i/ 2 / 4 (2) V 2 0( 2 ) V 2 ^ t (p . ks| 



2 2 4 3a 2 H 2>a 2 H 2 

Using Eqs. (6.8) of Ref. 0, together with Eqs. (|B8|I given in the present section, one obtains for the time evolution 
of the intrinsic second order terms <p^ 2 ' and tp^ the result 

V ' ; 3a 2 # 3a 2 # 9a 2 tf 9 a 4 H 3 K ' 

Now this is the first check we wanted to do. We can compare with Eq. I|39|l and see that deriving by parts and using 
aH = 2/t, they coincide (up to four-derivative total spatial gradients). 

An important observation here is that the presence of the intrinsically second-order terms is necessary to obtain not 
only the correct result, but even a consistent calculation. In fact if one ignores them, i.e., ignores ij}^ an d X m the 
synchronous gauge and ignores and ip^ in the Poisson gauge, and compares the two results, one sees that they 



Note however that what is "intrinsic" second order in a gauge, can become first order times first order in another one. 
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do not coincide. In other words they are not really "intrinsic" second order, since what appears as intrinsic second 
order in one gauge, can appear as first order times first order in another gauge. 

At this point we may also average expression Eq. (|B19J| over a volume d 3 x as we did in Sec. IIIII This has to be 
done with the appropriate measure of integration, which is d 3 Xy^y, where 7 = det(gij). Up to first order the spatial 
metric is given by 



dy a dyP . 



dx l dxi 



gap = dp 



3a 2 H 2 



did j if 



8ij{l-2<p) 



3a 2 H 2 



(B20) 



where g is the metric in the z, y coordinates, and where we used the transformation of Eq. (|B16|) . Using this metric, 
the measure becomes 



d 3 yy/j = d 3 ya 3 



l-5tp- 



3a 2 H 2 



(B21) 



and this is the same as Eq. (jfJJ), and so the final result is the same. 

Finally, we can check the computation of Ref. which used only the first order metric in this gauge. That is, we 
have to put (j)^ and ip^ to zero in Eq. (|B18|) . 



6{t,x l ) = 3H- 
= 3H- 



4 LD>3 k 



V <P,j k 



2VV 8yV 2 if 

3a 2 H ~ 9a 2 H ' A5a 2 H 9 a 4 H 3 ' 
2VV 128^>V,i 89i(9j9V) ^V' jk ^,jk 



3a 2 H 45a 2 H 9a 2 H 9 a 4 H 3 

and average it with the measure of Eq. i|B21|) . Doing so, we recover the result 

22{9 i <pd i (p) , 22(9i(^V)) , 



(9)(t)=3H- 



45 a 2 H 



9 



a 2 H 



(B22) 



(B23) 



as in Eq. (11) of Ref. 0. 
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